Non-Gaussianity-induced Heat Conduction 
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We study the properties of a Brownian motor that is driven by non-Gaussian noises from athermal 
environments. We find that new terms should be added to the conventional Fourier law and the 
fluctuation theorem for the heat current, where its average and fluctuation are determined not only 
by the noise intensities but also by the non-Gaussian nature of the noises. Our result explicitly 
shows the absence of the zeroth law of thermodynamics in athermal systems. 
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Introduction.- Recent experimental developments in 
single-molecule manipulation have enabled us to inves- 
tigate the detailed thermodynamic properties of small 
systems such as colloidal and biological systems [l[ . This 
direction of researches has shed new light on Brown- 
ian ratchets [I-14|, where the average and the fluctu- 
ation of the heat current are respectively characterized 
by the Fourier law [15|, [l6| and the heat fluctuation theo- 
rem 



17 



- 26[ . If the environments are in thermal equilib- 
rium, stochastic thermodynamics with Gaussian noises 
is a powerful methodology 27]-3(| to investigate non- 
equilibrium thermodynamics of Brownian ratchets. 

On the other hand, the effects of non-Gaussian noises 
from athermal environments have been reported in, for 
example, electrical circuits 3l| and biomolecular sys- 
tems [321 



The conventional approaches in stochastic 
thermodynamics are not straightforwardly applicable to 
such systems, because the environments are not in ther- 
mal equilibrium. One approach to this problem has been 
the use of non-Gaussian stochastic models 
We then raise the following question: how should the fun- 
damental thermodynamic relations, such as the Fourier 
law and the fluctuation theorem, be modified in the pres- 
ence of non-Gaussian noises? 

Our aim in this Letter is to answer this question. On 
the basis of non-Gaussian stochastic energetics [Ifjj], we 
study a Brownian motor driven by non-Gaussian noises 
from two athermal environments, and derive a general- 
ized Fourier law. The average heat current between the 
environments is determined not only by the difference in 
the noise intensities (i.e., the temperatures for the case 
of equilibrium environments) , but also by the difference 
in the non-Gaussianity of the noises. In particular, even 
when the noise intensities of the environments are the 
same, the heat can be conducted purely by the effect of 
the non-Gaussianity. We also derive a correction to the 
heat fluctuation theorem, which reveals the fundamental 
properties of the heat fluctuations in the presence of the 
non-Gaussian noises. Moreover, we investigate the valid- 
ity of the zeroth law of thermodynamics for athermal sys- 
tems, and find that the zeroth law is not universally valid 




FIG. 1: A schematic picture of the non-equilibrium Brownian 
motor. The noise from the left (right) environment is char- 
acterized by the noise intensity T (T') and the fourth-order 
cumulant K (K'). K and K' characterize the non-Gaussianity 
of the noises. 



but depends on the details of a contact device between 
two systems. We numerically verify our result, which 
demonstrates that the direction of the average heat cur- 
rent depends on the non-Gaussianity of the noises, and 
that the properties of the heat fluctuation significantly 
deviates from those of the conventional heat fluctuation 
theorem. Our result would serve as a theoretical foun- 
dation to study the energy transport and the irreversible 
phenomena in athermal systems such as electrical circuits 
and biological systems. 

Model. - We consider a non-Gaussian stochastic model 
of a Brownian motor which consists of two vanes that are 
attached to two athermal environments by a spring (see 
Fig. 1). The vanes are driven by athermal fluctuations 
in the environments, and the spring conducts an energy 
current induced by the fluctuations. We refer to the en- 
ergy current as the heat current. The motion of the vanes 
is described by the following Langevin equations: 



dx dU .„ „. - 

dy dU 
Tt=-^ v) + 7h 



(1) 
(2) 



where x,y are the angles of the vanes, U(x,y) is the 
dimensionless potential energy of the spring, and £, fj 
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are independent white non-Gaussian noises that char- 
acterize the fluctuations from the athermal environ- 
ments. In the following, (A) denotes the ensemble av- 
erage of a stochastic variable A, and the Boltzmann con- 
stant is taken to be unity. The first, second, and the 
fourth cumulants of the noises are given by (£(£)} = 

W)) = 0, = 2T5 (*i - *a), (v{h)v{t2)) = 

2T'6(t 1 -t 2 ), (i(ti)i(t 2 )i(t 3 )i(U)) c = K8 4 {t t ,t 2 ,hM), 
and {^{ti)f](t 2 )f]{h)f]{ti)) c = K'S i {t l7 t 27 t 37 t 4: ), where 
(£(ti) • ■ • £(t n )) c denotes a n-th order cumulant and 
5 n (t\, . . . , t n ) is a n-point delta function with an positive 
integer n [4J|. In the following arguments, the cumulants 
(£(ti) • ■ • t;{t n ))c with n = 3 and n > 5 do not contribute 
to our result. On the basis of stochastic energetics [27 
I30L l36| | , the heat current is defined by 



dQ ( dx A 



(3) 



where * describes a stochastic integral that is defined 
as a white noise limit of a colored noise. As discussed in 
Ref. [361 ] in detail, the * integral for an arbitrary function 
f(x(s)) is given by 



ds£(s) * f(x(s)) 



N-l 

lim lim N 

e->+0 At->+0 

i=0 



(4) 

where At = t/iV, tj = iAt, and £ e (ij) is the colored noise 
with a correlation time e > such that lim e _>o £e(£i) — 
£(ti). For non-Gaussian noises, the definition of the heat 
is not consistent with the Stratonovich integral but with 
the * integral (3(|. We note that the Stratonovich and * 
integrals are the same for Gaussian noises. In the case of 
K = K' = 0. we can show that the motor obeys the con- 
ventional Fourier law and the heat fluctuation theorem: 



J = kAT, 



lim l ln P(Q(t) = +9 *) 

t^oo t p(Q( t ) = -qt) 



(5) 



(6) 



where J = (dQ/dt)ss = ^^it^oo{dQ/dt) is the average 
heat current in the steady state, k is the thermal con- 
ductivity, AT = T - T', A/3 = 1/T - 1/T', and q is the 
time average heat current. 

Generalized Fourier's law. - We now discuss the gener- 
alized Fourier law in the case of a non-harmonic spring 
potential U(x,y) — {l/2)(x~y) 2 + (e/A)(x~y) A , where e 
is a small positive number. By using the first order per- 
turbation in terms of e, we obtain the generalized Fourier 
law: 



J = kAT + k'AK, 



(7) 



where k = 
AK = K 



(l/2)[l + {3e(T + T')/2}], k' = e/8, and 
- K' . This is the first main result of this 
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FIG. 2: (Color online) The numerical verification of Eq. ([7|. 
The dashed line is theoretically obtained from Eq. (JT)), the 
cross points show the numerical data of our Monte Carlo 
simulation, and the open circle is the point at which the 
heat current is switched. As the non-Gaussian property 
1/A becomes larger, the heat current becomes smaller. We 
assume the ergodic property for the system {dQ/dt)ss = 
lmiT-).oo \^-/T Jq ds (dQ/dtJ^ and calculated the long time 
average instead of the ensemble average. The time step is 
given by 0.001 and the entire time interval for the average is 
5 x 10 7 . The limit 1/A — > +0 is the Gaussian limit. 



Letter, which will be shown later. The first term on 
the right hand side (rhs) of Eq. ([7]) corresponds to the 
conventional Fourier law, and the second term describes 
the correction term due to the non-Gaussianity of the 
noises. This result implies that the heat is conducted 
from the environment with the higher non-Gaussianity 
to the other environment. We note that, if the potential 
is harmonic with e = 0, the effect of the non-Gaussianity 
vanishes in Eq. ([7|). 

We numerically demonstrate the validity of Eq. (JT]). 
For simplicity, we assume that is a white Gaussian 
noise and that f](t) is a two-sided Poisson noise with in- 
tensities ±7 and a transition rate A/2: 



V 



(8) 



where ti,§i are times at which Poisson flights happen. 
Figure 2 shows our numerical results with T = 0.30, 
AT = 0.01, and e = 0.1. We plot the average heat cur- 
rent by changing 1/A. The direction of the heat current is 
changed at 1/A* = (AT/eT 2 ) [1 + 3e (T + T) 2] ~ 1.24. 

We focus on the case of T = T, where AT on the 
rhs of Eq. ([7]) vanishes but AK drives the heat current. 
From 2T = I' 2 X, K' = 4T 2 /A and K = 0, we obtain the 
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average heat current 



J = - 



~2\ 



(9) 



where the rhs describes the effect of the non-Gaussianity. 
In the limit 1/A — > +0, the non-Gaussian contribu- 
tion vanishes, because the Poisson noise converges to the 
Gaussian noise due to the central limit theorem. 

Correction to the heat fluctuation theorem.- We next 
discuss a correction term to the conventional heat fluctu- 
ation theorem on the basis of the perturbative expansion 
in terms of q, A/3, K and K 1 . We assume that these four 
quantities are proportional to a small positive parameter 
6. To make the problem simpler, we consider the case 
of a harmonic potential with e = 0. We obtain a correc- 
tion term to the heat fluctuation theorem up to the third 
order perturbation in terms of 5: 



lim 



t P{Q(t) = 



-qt) 



-qt) 



q&P 1 



K + K' 
16T 2 



(10) 



which will be shown later. This is the second main result 
of this Letter. Although the conventional Fourier law ([5]) 
holds for a harmonic potential with e = 0, the conven- 
tional heat fluctuation theorem ([6]) should be modified as 
Eq. (|10[) even for the harmonic potential, which implies 
that the effect of the non-Gaussianity appears only in the 
higher cumulants in this case. 

We also numerically checked the validity of Eq. (fTU|) by 
using the two-sided Poisson model ©, as shown in Fig. 
3. By taking t = 1000, 1/A = 0.5, T = 0.20, and AT = 
0.01, we numerically obtain the fluctuation function 
/(<?) EE lim t ^ 00 (l/t)lnP(Q(t) = +qt)/P(Q(t) = -qt) 
and compare it with Eqs. (|10l) and ©. We observe a 
significant deviation from the conventional heat fluctu- 
ation theorem ([6]), and the deviation is consistent with 
our result (fT0|) . 

The zeroth law of thermodynamics. - We next discuss 
the zeroth law of thermodynamics 
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In the case 

of non-Gaussian noises with AK ^ 0, Eq. ([7]) implies 
that the condition of J = depends on the potential 
shape of the spring characterized by e. In contrast, 
in the case of Gaussian noises, J = if and only if 
AT = 0. Therefore, the zeroth law of thermodynam- 
ics is not universally valid for non-Gaussian noises; the 
condition of J = depends on the details of the contact 
device (i.e., the spring). When we fix the contact de- 
vice, however, the transitivity of equilibration holds and 
we can introduce an indicator characterizing the direc- 
tion of heat current. We note that a similar observa- 
tion was reported for driven lattice gas systems iol 41|. 
To show this, we consider three athermal environments, 
AEi, AE2 and AE3, whose fluctuations are characterized 
by the second and forth order cumulants (Ti,Ki) with 
i = 1,2,3. If we link the contact device between AE^ and 
AEj (i,j = 1,2, 3, i 7^ j), the average heat current be- 



tween them is given by 



*>ij A^jT 
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FIG. 3: (Color online) The numerical verification of Eq. (|10[) . 
The red broken line is obtained from the generalized fluc- 
tuation theorem (|10|l . the black chain line is obtained from 
the conventional fluctuation theorem ©, and the blue cross 
points show the numerical data of our simulation. We adopted 
the Monte Carlo method to make the histogram of the heat 
distribution function, and obtain the numerical data of the 
fluctuation function f(q). The bin- width for the heat his- 
togram is 0.0001, the time step is 0.0002, and the number of 
samples in the Monte Carlo is 5 x 10 7 . Due to the large cost 
for the calculation, the numerical data of f(q) is restricted to 
|g|<0.2. 



Kij ee (l/2)[l+{3e(T+T J )/2}], ee 



AijT EE T; 



and AijK 



We can then easily show the tran- 
sitivity of equilibration: if J12 = and J23 = 0, then 
J13 = 0. We also introduce a device-dependent indicator 
fijT, K) as 



1. 



3e, 



fr(T,K) = -T+—T 



-K, 



(11) 



which characterizes the direction of the average heat 
current = /^(Tj, .£Q) — [i e (Tj, Kj). In this sense, 
H e (Ti, Ki) plays a similar role to that of the temperature 
in equilibrium thermodynamics. 

Derivations of the main results.- We now prove 
Eq. (J7]). By introducing a new variable z = x — y, Eqs. ([1]) 
and ((2]) reduce to a single equation: 



— = -22 - 2ez 3 +£-f). 



For small 
expanded 

z {t) = 

-2j>e-*C« 
culation 



as 
t 



dz 
~dt 

the 



(12) 



solution 

ZQ + 

s) z 3 (s). By 
we obtain (z 2 )ss 



of Eq. {T5J can be 
ez\ + 0(e 2 ), where 

-fj(sfj and zi(t) = 

a straightforward cal- 
= (T + T')/2 - 



where 



3e(T + T') 2 /4 - e(K + K')/8 + 0(e 2 ), (z 4 ) SS = 
3(T + T') 2 /4+(K + K')/8 + 0(e), (5*C)ss = T + 0(e 2 ), 
and (z 3 * |) ss = 3T(T + T')/2 + if/4 + O(e). The 
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corresponding stochastic heat current is given by 



dQ/dt 
obtain 

J = 



2ez 4 



z * £ + ez A * £ + 0(e 2 ). We then 



-<z 2 ) 



ss 



2e(£ 4 )ss + (^£)ss + e(^*£)ss (13) 



1 



3e(T + T') 



AT 



-AK, 



(14) 



which implies Eq. ([7]). 

To derive the generalized heat fluctuation theo- 
rem (ITU1) . we set e = and consider the large deviation 
principle liHI^ 



lim -In P(Q(t) = qt) = v*q+ (15) 

t— S- + 0O £ 

where Vt'(w) = lim t ^ +tx) In (e~ v ®W) /t is the cumulant 
generating function, and v* is the point at which vq + 
^(v) takes the minimum. For the Gaussian noise, the 
cumulant generating function is exactly given by [431 ] 



*o(») = 1 - y/(l-Tv)(l + T'v). 



(16) 



For the non-Gaussian noises, the cumulant generat- 
ing function can be expanded as &(v) — ^o(v) + 
&i(v) + 0(S 2 ), where ^i(v) is the order of 5. By 
a straightforward calculation, we obtain the variance 
of the heat flux a 2 (t) = ((Q(t) - (Q(t))) 2 )/f = 
(T + T'f/At + (K + K')/m + 0{t- 2 ), where Q(t) = 

J * ds (^— z 2 (s) + z(s) * £(s)J and z(t) = z a (t) in the case 
of e = 0. We then obtain 



' ,,2 



K + K'v 
16" 



0{Sv 3 



(17) 



where we consider up to the second order term of v be- 
cause v* is the order of 5. From Eqs. (fT5)) . (fTo) . and (TT71) . 
we obtain 



1 



lim -\nP(Q(t)=qt) = - 
t-xx t 21 1 



q 2 + qAT + AT 2 1 A 



{K + K')/8 



0(5% 
(18) 

which implies the generalized fluctuation theorem (I10[) . 

Conclusion and discussion.- In this Letter, we have 
studied a Brownian motor driven by non-Gaussian noises 
from two athermal environments. As a result, we found 
new additional terms in the Fourier law and the heat fluc- 
tuation theorem, which implies that the heat current can 
be induced by the non-Gaussianity of athermal fluctua- 
tions. We also discussed that the zeroth law of thermody- 
namics is not universally valid for athermal systems. Our 
numerical results are not consistent with the conventional 
Fourier law and the fluctuation theorem, but consistent 
with our analytical results obtained in this Letter. Our 
results would open up an extension of nonequilibrium 
thermodynamics to athermal systems with non-Gaussian 
noises. Our theory can be experimentally tested in elec- 
tric circuits and biological systems driven by the Poisson 
noises [3ll. |32|. 
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SUPPLEMENTARY INFORMATION: 
CUMULANT FUNCTIONAL AND n-POINT 
DELTA FUNCTIONS 

Here we discuss the relationship between a cumulant 
functional of the non-Gaussian noise £(t) and a n-points 
delta function. A cumulant functional and a n-th 
order cumulant (£(ti) . . . £(i„)) c [45j are respectively in- 
troduced as 



<£(tl)--.£(tn)> 



$H ee In (e^o^M*)^ 
6 n <S>\v] 



5iv(t\) . . . Siv(t n ) 



(19) 
(20) 



v=0 



where v(s) is an arbitrary function. According to the Ito- 
Levy decomposition 4(| 47 1, the cumulant functional can 
be transformed into the standard form of Levy processes 



ds 



iav(s) - 



a 2 v 2 {s) 



dz 



(21) 

where a and a 2 are arbitrary constants, w(z) is a transi- 
tion rate function. We next introduce the n-points delta 
function as 



5 n (ti, . . . , t n ) 



oo (h = ■■■ = t n ) 



(else) 
dt 2 . . .dt n S n (ti, ...t n ) = l. 



(22) 
(23) 



We assume the symmetric property for the delta func- 
tion S„(ti,...t n ) = 5 n (si,...,s n ), where {si,...s„} is 
an arbitrary permutation of {ti, . . . , t n }. The following 
equality is crucial to derive Eqs. (7) and (10) of the main 
text: 



dU 



(t, t2, ■ ■ ■ , t n ) 



1 

n 



(24) 



The derivation of this equality is as follows. By definition, 
we obtain 



dti . . . dt n S n (ti, . . . , t n ) — t. 



(25) 
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By differentiating the left and right hand side around t, 
we obtain 



dt 2 ...dt n 5 n (t,t 2 ,t 3 ,...,t n )+ (26) 



+ / dhdt 3 ...dt n S n (t u t,t 3 ,...,t n ) + ... (27) 



The symmetric property of the delta function leads to 
the following equality 



n dt 2 ■ ■ .dt n S n (t,t 2 , ■ ■ ■ ,t n ) = 1, 



(29) 



which implies Eq. (|24[). We can represent the functional 
derivatives of the cumulant functional with the n-points 
delta functions: 



Siv(ti)Siv(t 2 ) 



v=0 



d 4 $[v] 



8iv{ti)8iv(t 2 )5iv{t 3 )8iv{ti) 



2T6 2 (t u t 2 ), (30) 



Ki8i(t\,t2, h, t 4 ), 



[12 

[is; 

[14 

[is; 

+ / dt 1 ...dt n ^ 1 6 n (t 1 ,t 2 ,...,t) = 1. (28) [16 



[ir 
[is; 

[19 
[20 

[21 
[22 
[23 

[24 
[25 

[26 

[27 
[28 
[29 
[30 

[31 

[32 
[33 

[34 

[35 

[36 

[37 

[38 
[39 
[40 

[41 

[42 

[43 
[44 
[45 
[46 



(31) 



where the second and fourth order cumulants are re- 
spectively introduced as 2T = a 2 + J z 2 w(z)dz and 
K = J z A w{z)dz. We note that the assumption of the 
symmetry of the delta function is consistent with that of 
the mixed functional derivative 



5 n <f>[v] 



Siv(t\) . . . Siv(t n ) 



5iv(si) . . . Siv(s n ) 



(32) 



v=Q 
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